Closed form solution for the self-resonances in a short 

Josephson junction 



b,* 



S. De Nicola% M. Adamo^, E. Sarnelli'^, C. Nappi 

""CNR Istituto Nazionale di Ottica 
1-80078, Pozzuoli, Napoli, Italy 
^ CNR Istituto di Cibernetica "E. Caianiello" 
1-80078, Pozzuoli, Napoli, Italy 



Abstract 

We present a closed form solution for the self-resonances in a short Josephson 
;unnel junction. This solution is alternative to the well known textbook result 
jj, |2|] based on a series expansion. Results are derived for the up-to-date case 

of a — TT junction. 
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When a constant voltage V is present across the electrodes of a Josephson 
junction, the current flowing into the junction oscillates at a frequency u = 
2'kV/^q, where $o = 2.07 x 10"^^ Weber is the flux quantum. On the 
other hand, if the major size of the junction L is shorter than the Josephson 
penetration depth Xj (short junction limit), the junction can be viewed as 
a cavity of length L. In the presence of a uniform magnetic field and at 
an applied junction voltage Vn = $oi^n/27r, the oscillation frequency of the 
Josephson current matches the n-th harmonic of the junction cavity mode 
resulting in the excitation of some of the modes at the frequencies Un/'iTc = 
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nc/2L (n = 1,2,3...), where c is the hght velocity of the waves in the resonator. 
In th. ca.e typ.cal cu.e,. .tep., known a. T.ie .tap. y, appea. .n the 
I — V characteristic of the junction pj. 
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Figure 1: Schematic of the superconductor- insulator- ferromagnetic metal-superconductor 
— TT Josephson tmmcl junction, along with the coordinate system used in this work 



Kuhk [2j] showed how to calculate, in agreement with the experiments 
jsj], the dependence of the amplitude of these steps as a function of an ap- 
plied magnetic field in the short junction limit. In the Kulik's solution, the 
Josephson phase difference is expressed in the form of a Fourier series, each 
term representing one of the resonances at the voltages Vn- Fiske steps are 
observed also in long junctions (L > Aj), when a magnetic field is externally 
applied. In this case however the junction behaves more like a transmis- 
sion line than a resonator and relevant to the underlaying mechanism is the 
presence of so called /Jiuxons, particle-like current-field structures periodically 
driven inside the junction 4|. 

In this letter we present and discuss an alternative equivalent approach 
we have developed for calculating the Fiske resonances in a short Josephson 



junction. The method is based on the development of a closed form solution 
to the linearized sine-Gordon equation. Compared to the standard result this 
method may present some advantages, in particular if one is interested to a 
description of the phase dynamics given in closed form. We apply the method 
to the relevant case of a — tt Josephson junction, i.e. a junction which has a 
coupling, changing between and vr along the junction length, which implies 
an anomalous current-phase relation in the vr- region (see equation([T])). This 
physical situation can be realized, for instance, in superconductor-insulator- 
ferromagnetic metal-superconductor tunnel junctions, i.e. junctions in which 
a ferromagnetic layer with step-like thickness is inserted, in addition to an 
insulating layer [5|, la] (see figure [1]) . 

The effect of the presence of few, or many, adjacent — tt phase shifts on 
the self resonant modes of a Josephson junction has been considered in the 
context of YBCO grain boundary Josephson junctions j?]. In that reference, 
in order to calculate the contribution of the self-resonances to the current- 
voltage characteristics, an extension of the Kulik expansion was developed. 
The results have been successfully used to fit data on Fiske steps in — vr 
Josephson tunnel junctions [8]. Also recently, a careful analysis of Fiske 
modes, based on the Kulik theory, has been carried out in superconductor- 
insulator-ferromagnetic metal-superconductor to extract information on the 
junction quality factor and the relevant damping mechanisms 9(]. 

Let us consider a — tt junction with two adjacent regions [0, /] and [/, L], 
characterized by two different maximum Josephson current densities, ji and 
j2 (ji > and j2 > 0) and a tt phase shift in region [I, L] (See figure [1]). The 
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supercurrent density can be written as 

{?i sin(o9) if < w < / 
(1) 
j2 sm{(p + Tc) if I < y < L 

Assuming a one- dimensional system {w << L), the equation for the phase 
difference f{y,t) is expressed by 



1 1 dip Jc{y) ■ ^ , nf ni Je 

W ~ ~ = 1^ + '"'^ - "I - (X> 

where we have introduced the specific normal resistance R and capacitance 
C of the junction. Function 9{y) is the Heaviside step function, je/{Jc) is 
the normalized external bias current density, and (Jc) = (jV + j2{L — I)) / L 
is the average maximum Josephson current density. In Eq. [2] the aver- 
age Josephson penetration depth is given by {Xj) = a/^'o/ (2vr/io(Jc)'^e//) 
7| where the effective depth is d^fj ~ 2Al with the London penetra- 
tion length and we have introduced the average plasma frequency (up) = 
V27r(Je)/$oC. 

In the geometry of figure [T|, the magnetic field within the junction is 
related to the derivative of if through the relationship 

Hiv)--.^('4) (3) 



27r/io4// \dy 

The boundary conditions appropriate to equation ([2]), and to the considered 
geometry are 



dip 
dy 



dip 
dy 



2TrflodeffHe 
$0 



4 



Further important conditions are given by the continuity of and of its 
derivative at the point y = I. The last condition expresses the continuity 
of the magnetic field. From now on, all lengths and time variables will 
be normalized to (Aj) and to the inverse average plasma frequency (wp), 
respectively, so that equations ^ and (jl]) become respectively 

^yy - ^tt - a^Pt = ^jjj sin [(p + 7re{y - /)] - (5) 



(fy{0,t) = 'fyiL,t) = -he, (6) 

where a = 1/Q = l/{(jjp)RC is the dimensionless damping coefficient, Q 
the damping quality factor and h^. = Hf./ {Xj){Jc) is the normalized external 
magnetic field. 

Following the Kulik approximation, for solving equation ([2]) we write the 
phase ip as a sum of two terms p{y,t) = ipQ{y,t) + {pi{y,t) where the un- 
perturbed term is po = ut ~ hey, pi ^ po and v = ^q/ 271 dipi/dt is the 
perturbation to the steady voltage V. Here u = {2'kV/ $o) / {^p) is the nor- 
malized Josephson frequency corresponding to the fixed voltage V applied 
between the electrodes of the junction. We note that, to the zero-th order, no 
magnetic field is associated with the 0— vr discontinuity in the present approx- 
imation, as on the left and right of the point y = I we have dpo/dy = —he, 
and the only magnetic field present is the external one. A time dependent 
magnetic field perturbation however appears to the 1st order, i.e. dpi/dy. 

If we denote the perturbation pi as p^, when considered in the interval 
< y < ^ and pf, when considered in the interval I < y < L, the linear 
equations providing p^ and pf are respectively 
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for < ?/ < /, and 

ioT I < y < L. In equations (I7j) and (|H]) we have defined gi = ji/{Jc), 
g2 = —j2/{Jc)- The boundary conditions appropriate to equations ([7]) and 
dHD are 

dy dy 

¥,f(/,t) = y,f(/,t) (10) 

dy dy 

Equation ([9]) is the requirement of perfect reflectivity of the edges of the 
junction and continuity of the phase and of its first derivative at the point 
y = I are determined by the equations ffTOl) and f|TT]) . respectively. 

After defining the two complex functions u{y) and v{y) through the fac- 
torization 

V?[ = Re [u{y) exp(-iwt)] (12) 

(^f = Re [v{y) exp{-iut)] (13) 

and substituting (fT2!) and (fT3|) in ([7j) and ([8]), we find that u and v satisfy 
the two equations 

— j^^^u = igi exp{ihey) (14) 



for < y < /, and 

— + ^^v = ig2 exp{ihey) (15) 

for I < y < L, where we have defined = (u'^ + iau). The boundary 
conditions for the two complex functions u{y) and v{y) are 

dy dy 

nil) = vil) (17) 

du{l) ^ dvjl) 
dy dy 

We can write now the general solution to equations ( !T4l) and ( !T5l) in the 
following form 

u{y) = a^e-'^y + /3„e+*^^ + Mp(i/), <y <l (19) 



v{y) = a,e-'^^y^ + /3,e+^>^(^) + Vj,{y),l <y<L (20) 

where the basic task reduces to determine the four constants /S^, ay , 

I3v from the boundary conditions (|T6 l) -( fT8l) . In equations ( fT9l) and ( l20l) we 

have introduced the particular solutions Mp(?/) and Vp{y) which, by following 
standard methods, can be promptly written as 

My) = (21) 

X - hi 
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By using the boundary conditions (I16l) - fll8l) we find the following relationship 
between the unknown coefficients a„, and (3v 

iX (a« - Pu) = 2^\o (23) 



x^-K 

IX (a« -a^)e^ + ix [Pu - Pv) e ^ = 5 

X - hi 

From equations (PT]) - (B51) we obtain au, f3u, cn^, which can be written as 

"'^ = 2x{smxL){x'-hl) ^^^^ 
[{92 - 9i)e"'^'S{L + Kig.e-'^"^ - g^e'^^"^^ 



2x(sinxL)(x2-/^2) ^ (2^) 



= TT^ T\( 2 h2^ ^ (27) 

2x(sin xL){x - hi) 

[{92 - gi)e'^^^'^^'^S{-l) + /ie(^7ie-*^^ - ^726^'^^ 
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where S{z) = he cos (x-z) + ix sin (x-^)- 

Equations and ([20]), with coefficients (IMj). (|25|). (|26|) and 

[27]), determine in the present approximation of 'small Q' and short junction 
l|, the dynamics of the phase ipi and the magnetic field inside the junction 

dipi/dy for arbitrary u and values. 

Next, in order to extract a possible dc term in the current, we have to 

carry out time and space averages. That is to say, we have to calculate the 

quantity 



Jdci'^, he 
-L 



1 



- / {Je{y) sin ip{y))dy 



(28) 



1 
L 
1 
L 



{Je{y) sm{ipo{y) + ipi{y, t) + n9{l - L)))dy 







Ji 



(sin(v9o + '^{))dy - j2 / (sin(v?o + '^i))dy 



Angle brackets indicate time average over the period T = 27r/a;, i.e. if 
F{t) is an arbitrary function of the time, then {F{t)) = (1/T) F{t)dt = 
(l/27r) J^"^ F{ujt)d{ut). Furthermore, since (sin((y9o+v5i)) ~ (siny9o) + (v^i cosy^o) 
and 

(sin(/?o) = (29) 
(<^fcos<^o) = ^i?e [ue-'^'^y] 

(</.fcos</.o) = ^i?e [ve-^'^^y] 

we obtain for the dc component of the current due to the self- resonances, 
the expression 
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CL>l«jJp> 

Figure 2: Current density vs uj in zero external magnetic field in a — tt Josephson 
junction; a — 1] L = l\he = 0]l — 0.5; g\ = 1; 52 = —1; self-resonances appear only at the 
odd positions w = tt, Stt, Stt 

The result can be expressed in term of the normahzed flux applied to the 
junction, with the position heL = 27r0 = 27r$/$0; where $ = fioHedeff{L 
and deff and L are expressed in the usual units in the last expression). 
The coefficients, equations (124|) . (126|) and (127|) . diverge for sin^L = in 
the limit of vanishing damping. This gives the resonance frequencies of the 
system Un = nn/L and, with very good approximation, the frequencies at 
which, in the presence of damping, the current equation ( l30l) peaks. The 
amplitude dependence of the n—th step on the magnetic field, of fundamental 
importance for a comparison with the experiments , can be calculated 
by equation f l5U]) by setting the value oi u aX u = Un- The same equation 
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Figure 3: Current density vs uj in the presence of an external magnetic field in a — tt 
Josephson junction. The two curves refer to two different values of the magnetic field 
(/le = f f and hf, — f 0, thin and thick line respectively); a = 1; L = 1; I ^ 0.5; gi = f ; g2 = 

-1; 

can be also used to probe the 'shape' of the resonances as a function of the 
frequency and Q damping factor. For the sake of illustration, in figures [2] and 
iniwe show two cases of current versus frequency obtained by using equation 
f pOj) . The first graph refers to the case of zero external field. As can be seen, 
odd resonances persist in zero field, even though the amplitudes of those 
following the first are vanishingly small, a phenomenon typical of the — vr 
junction ^|. The second one refers to a generic situation of presence of an 
external magnetic field and the two curves are calculated at two different 
values of the normalized field. The right hand side of each of the bell shaped 
peak has negative resistance and, for this reason, has no relevance for a 
comparison with the experiments, where usually a current bias set up is 
considered. We point out that, in principle, in the framework of the Kulik 
theory to obtain the same accuracy in the determination of the current 
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density as a function of frequency or magnetic field, one would have to sum 
up the contributions of the entire series representing the current. Finally it is 
worthwhile to stress that the result for a — tunnel junction, with uniform 
maximum Josephson current density, can be easily recovered by the above 
method. In this case we have to discuss only the equation 

+ X^uo = i exp{ihey) (31) 

ay 

with < y < L, ipi = Re[uo{y)exp{—iujt)] and the boundary conditions 
duQ{0)/dy = duQ{L)/dy = 0. The solution can be written in the following 
form 

uoiy) = aoe-'^y + + (32) 

X '^e 

only two coefficients ao e /3o have to be determined. These coefficients are 
obtained from the two boundary conditions since, now, the continuity condi- 
tion at / = is no longer required and they are given by the following simple 
expressions 

"° 2x{sui xL) {x^ - hi) ^ ^ ^^"^^ 

2xismxL)ix'-hl)^ ^ 
It is easy to verify that these coefficients can be formally obtained by ( 12^ - 
(!27|) by letting / = and taking gi = g2 = ^- In the same limit, we can also 
determine from the general result given by equation (1301) . the explicit expres- 
sion of the frequency and field dependance of the current density, namely 



Jdc{u}, he) = Re 



ihl 

(cos xL — cos h(,L) + 



xL{smxL){x^ - hl)"^ 2{x^ - h^) 



(34) 
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In conclusion, we have presented and discussed a closed form solution for 
the determination of the dynamics of the phase in a Josephson junction in 
the limit of short junction. Within this framework we have also derived an 
expression for the dc current density associated to the Fiske resonances in a 
0— TT Josephson tunnel junction. This approach can be relevant for improving 
the accuracy of data fitting in the determination of the damping mechanism 
in — TT Josephson junctions or conventional Josephson tunnel junctions. 
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